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ABSTRACT 

This paper provides the first study of a new dynamical instability in superflu¬ 
ids. This instability is similar to the two-stream instability known to operate 
in plasmas. It is analogous to the Kelvin-Helmholtz instability, but has the dis¬ 
tinguishing feature that the two fluids are interpenetrating. The instability sets 
in once the relative flow between the two components of the system reaches a 
critical level. Our analysis is based on the two-fluid equations that have been 
used to model the dynamics of the outer core of a neutron star, where superfluid 
neutrons are expected to coexist with superconducting protons and relativistic 
electrons. These equations are analogous to the standard Landau model for su¬ 
perfluid Helium. We study this instability for two different model problems. First 
we analyze a local dispersion relation for waves in a system where one fluid is 
at rest while the other flows at a constant rate. This provides a proof of prin¬ 
ciple of the existence of the two-stream instability for superfluids. Our second 
model problem concerns two rotating fluids confined within an infinitesimally 
thin spherical shell. The aim of this model is to assess whether the two-stream 
instability may be relevant (perhaps as a trigger mechanism) for pulsar glitches. 
Our results for this problem show that the entrainment effect could provide a 
sufficiently strong coupling for the instability to set in at a relative flow small 
enough to be astrophysically plausible. 


1 INTRODUCTION 


In this paper we describe a n ew dynamical instabi lity in superfluids. This “two-stream” instability is analogous to the 
Kelvin-Helmholtz instability ( Drazin fc Reid 1981 ). Its key distinguishing feature is that the two fluids are interpene¬ 
trating rather than in contact across an interface as in the standard scenario. The two-str eam instability is well known 
in plasma physics [where it is sometimes referred to as the “Farley-Buneman” instability (Farley 1963; Buneman [l963 


1959)1, and it has also been discussed in various astrophysical contexts like merging galaxies (Lovelace et aj 1997) and 


pulsar magnetopheres (|Cheng fc Ruderman 1977| ; Weatherall 1994; Lyubarsky 2002), but as far as we are aware it has 
not been previously considered for superfluids. In fact, the “standard” Kelvin-Helmholtz instability was only recently 
discussed in the context of superfluids ( |31aauwgeers et a] 2001). 

The similarity of the equations used in plasma physics [a nice pedagogical description of the plasma two-stream 
instability can be found in Anderson et aj (2001)] to the ones for two-fluid superfluid models inspired us to ask whether 
an analogous instability could be relevant for superfluids. That this ought to be the case seemed inevitable. To prove 
the veracity of this expectation, we have adapted the arguments from the plasma problem to the superfluid case, and 
discuss various aspects of the two-stream instability in this paper. 

Of particular interest to us is the possibility that the two-stream instability may operate in rotating superfluid 
neutron stars. Mature neutron stars are expected to be sufficiently cold (eg. below 10® K) that their interiors contain 
several superfluid/superconducting components. Such loosely coupled components are usually invoked to explain the 
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enigmatic pulsar glitches, sudden spin-up events where the observed spin rate jumps by as much as one part in 10® 
(Lyne et al 200C). Theoretical models for glitches (Eluderman 196!:; Baym et al 1969; Anderson & Itoh 1975) have 


been discussed ever since the first Vela pulsar glitch was observed in 1969 (Radakrishnan & Manchester 1969; Reichley 
& Downs 1969), but these event are still not well understood. After three decades of theoretical effort it is generally 
accepted that the glitches arise because a superfluid component can rotate at a rate different from that of the bulk of 
the star, and that a transfer of angular momentum from the superfluid to the crust of the star could lead to the observed 


phenomenon. The relaxation following the glitch is well explained in terms of vortex creep [see for example Cheng et 


al (198E)], but the mechanism that triggers the glitch event remains elusive. In this context, it seems plausible that the 
superfluid two-stream instability may turn out to be relevant. 


2 PROOF OF PRINCIPLE: A LOCAL ANALYSIS 
2.1 Superfluid dispersion relation 


We ta ke as our starting point the two-fluid equations that have been used to discuss the dynamics of superfluid neutron 
stars ( Mendell 1991a ; Andersson fc Comer 2001 ; Prix 2002| ). Hence, we consider superfluid neutrons (index n) coexisting 
with a conglomerate of charged components (index p). The constituents of the latter (mainly ions and electrons in the 
neutron star crust and protons and electrons in the core) are expected to be coupled by viscosity and the magnetic 
field on a very short timescale. Hence, we assume that these charged components will move together and that it is 
appropriate to treat them as a single fluid. 


The corresponding equations are (Andersson & Comer 
dtux + V • (nxvx) = 0 , 


2001 


Prix 


2002) 


( 1 ) 


where nx represent the respective number densities and vx are the two velocities. Here, and in the following, we use 
the constituent index X which can be either n or p. The respective mass densities are obviously given by px = rnx nx 
and we further introduce the relative velocity A between the two fluids as 

A = Vp-Vn ■ (2) 

The first law of thermodynamics is defined by the differential of the energy density or “equation of state”, £ = £!(pn, Pp, A^), 
namely 


d£ = p" dpn -I- p£ dpp + a dA^ , (3) 

which defines the chemical potentials and the “entrainment” function a as the thermodynamic conjugates to the 
densities and the relative velocity. With these definitions we can write the two Euler-type equations: 


(9t + Fn • V) (un -I-EiiA) -1- V (^ -1- p“) + enAjVw^ = 0, 

(4) 

{dt + Vp-V){yp — £pA^+V(^ + 'fP)—epAjVvl = 0, 

(5) 

where we have introduced the dimensionless entrainment parameters 


2a 

(6) 

ex = — . 

px 

The equation for the gravitational potential >1> is 


V^-l> = IttGp , 

(7) 


where p = pn -I- Pp. 

When a A ^ these equations make manifest the so-called entrainment effect. The entrainment arises because the 
bare neutrons (or protons) are ’’dressed” by a polarization cloud of nucleons comprised of both neutrons and protons. 
Since both types of nucleon contribute to the cloud the momentum of the neutrons, say, is modified so that it is a linear 
combination of the neutron and proton particle number density currents (the same is true for the proton momentum). 
This means that when one of the fluids begins to flow it will, through entrainment, induce a momentum in the other. 
Because of entrainment a portion of the protons (and electrons) will be pulled along with the superfluid neutrons that 
surround the vortices by means of which the superfluid mimics bulk rotation. This motion leads to magnetic fields being 
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The superfluid two-stream instability and pulsar glitches 3 


attached to the vortices, and dissipative scattering of electrons off of these magnetic fields. This “mutual friction” is 
expected to provide one of the main dissipative mechanisms in superfluid neutron star cores (Mendell 1991b|) . 

In order to establish the existence of the superfluid two-stream instability we consider the following model problem: 
Let the unperturbed configuration be such that the “protons” remain at rest, while the neutrons flow with a constant 
velocity vq. For simplicity, we neglect the coupling through entrainment , i.e. we take a = 0, and we also neglect 
perturbations in the gravitational potential. Under these assumptions, the two fluids are only coupled “chemically” 
through the equation of state. 

Writing the two velocities as Un = [uo + Svn{t,x)]x and Vp = Svp{t,x)x where Svn and (5up are taken to be suitably 
small, we get the perturbation equations 


dt3nn + vodxSuu + n^d^Svn = 0 , 

dt&rip -\- HpdxSvp = 0 , 

and 

dtSvn vodxSv^ -I- dxS'JT = 0 , 

dtSvp -I- dxShf = 0 . 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 


Next, we assume harmonic dependence on both t and x, i.e. we use the Fourier decomposition Svx{t, x) = 5vx exp[i(tjt — 
kx)] etcetera. This leads to the four equations 


i{uj — kvo)5n-n — ikn-nSv^ = 0 , 

iujSup — ikUpSvp = 0 , 

i{iM — kvo)Svn — ikS]T = 0 , 

iuj5vp — ik5yF = 0 . 


( 12 ) 

(13) 

(14) 

(15) 


We thus have four equations relating the six unknown variables Svx, Snx and Sfj,^. To close the system we need 
to provide an equation of state. Given an energy functional £ = £(nn,np) we have 




Kdrin ) 

np Xdnp) 


, d^£ , d^£ , 

^5 

(jTl-pCjTL'n. 


and similarly 

,~p d'^E , d'^E , 

mow ~ ~ ;; + TT-^OUp 


dupdrip 


dnl 


Finally, we define the two sound speeds by, cf. Andersson fc Comer (2001), 


dh" 

‘^pp 

drip 

m det S 


Tip 

dhF 

^p ‘^nn 

dup 

m det S 

Tip 


and introduce the “coupling parameter” 


C — Tin 


dJT 

drip 


drip 


Tin 


m det S 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


which also has the dimension of a veloc ity square d . For later convenience we have given the relation to the coefficients 
of the “structure matrix” Sxy used by Prix et a] ( 2002 ). 

With these definitions we get 

nnSjT = CpSrin -f CSrip 


UpS^jF = —CSrin + CpSrij 


p 1 


and we can rewrite our set of equations as 


rinSvn = 


~ ^ ~ florin+CSnp^ , 


( 21 ) 

( 22 ) 


(23) 
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ripSvp = = — 

k 


UJ 


—CSun + Cp(5nr 


Reshuffling we get 
2 


(i 

(i)‘ 


- 


(5nn = C(5np 


(5np = 


and a dispersion relation 


’//, \ 2 

/ ^ \ 2 


fUjC 2 

1 

1 

1 




-CSUn 




Introducing the “pattern speed” (the phase velocity) ap = uj/k we have 

\r 21 r 2 21 

[(cTp - Wo) - CnJ [(Jp - CpJ = - C . 


(24) 

(25) 

(26) 

(27) 

(28) 


Not surprisingly, this local dispersion relation is qualitatively similar to the one for the plasma problem (Anderson 
et ^ 2001). We will now use it to investigate under what circumstances we can have complex roots, i.e. a dynamical 
instability. 


2.2 The superfluid two-stream instability 

First of all, it is easy to see that (|^) leads to the simple roots 

_ J ±Cp 
^ ^ Wo ± Cn 


(29) 


in the uncoupled case, when C = 0. This establishes the interpretation of cx as the sound speeds. 

To investigate the coupled case, we introduce new variables x = (Jp/cn and y = vo/c-a- Then we get 

f{x,y) = ^[{x-yf - l][x^ - 6 ^] = 1 


where we have defined 

2 Thp 


Un Cn 


and 



(30) 


(31) 


The onset of dynamical instability typically corresponds to the merger of two real-frequency modes. If this is the 
case, a marginally stable configuration will be such that ( ^ ) has a double root. This happens when an inflexion point of 
f{x, y) coincides with /(r, y) = 1. This is a useful criterion for searching for the marginally stable modes of our system. 

As a “proof of principle” we consider the particular case of = 0.0249 and — 0.0379 (we will motivate this 
particular choice in Section HE). The real and imaginary parts of the mode-frequencies for these parameter values are 
shown as functions of y in Figure]^ We have complex roots (an instability) in the range 0.6 < y < 1.5. The corresponding 
mode frequencies lie in the range 0.03 < xo < 0.36. The fastest growing instability occurs for y « 1.1 for which we find 
that Im X « 0.17. In other words, in this particular case we encounter the two-stream instability once the rate of the 
background flow is increased beyond 


Wo = Cny « 0 . 6 cn 


(32) 


The corresponding frequency is given by 

UJ = kCnXo « O.lfcCn . (33) 


From this we see that the instability is present well before the neutron flow becomes “supersonic”. This is crucial since 
one would expect the superfluidity to be destroyed for supersonic flows. 

We have thus established that the two-stream instability may, in principle, operate in superfluids. Our example 
indicates the existence of a lower limit of the background flow for the instability. This turns out to be a generic feature. 
In contrast, the plasma two-stream instability can operate at arbitrarily slow flows. An ideal plasma is unstable to 
sufficiently long wavelengths for any given wq. In reality, however, one must also account for dissipative mechanisms. In 
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The superfluid two-stream instability and pulsar glitches 5 




Figure 1. Real (left panel) and imaginary (right panel) parts of the four roots of the dispersion relation (js^ for model parameters 
a? = 0.0249 and 6^ = 0.0379. For these parameters the quartic dispersion relation has four real roots for both y = 0 and y = 2, 
while it has two real roots and a complex conjugate pair for y in the range 0.6 < y < 1.5. In this range, the two-stream instability 
is operating. The grey area corresponds to 0 < Re a; < 1/2 which is contained in the range of the instability criteria discussed in 
Section IIC. 


the case of real plasmas one finds that the so-called Landau damping stabilizes the longest wavelengths (Anderson et 
al|001). Thus the two-stream instability sets in below a critical wavelength in more realistic plasma models, and there 


is typically (just like in the present case) a range of flows for which the instability is present. We will discuss the effects 
of dissipation on the superflnid two-stream instability briefly in Section IV. 


2.3 Necessary criteria for instability 


It is useful to consider whether we can derive a necessary condition for the two-stream instability.To approach this 
problem in fnll generality would likely be quite complicated, but we can make good progress for the simple one¬ 
dimensional toy problem discussed above. 

We begin by multiplying the Euler equation ( [l4| ) for the neutrons by the complex conjugate Sv^. This leads to 
(after also using the continuity equations to replace the perturbed number densities) 


ui — kvo — 


k Cn 


LO — kvo 


\&Vn\^ = C^-^Sv*Svp 

riuU) 


Similarly, we obtain from the second Euler equation di) 


7^2 ^2 
k Cn 


LO — ■ 


LO 


\ 5 vp\^ =C 


- - — 8vl5vn 

uj — kvQ 


Next we combine these two equations to get 
^2 


L = - (7- 


(Tp Vq 


( 7 p Vq 


Idunl^ + (o-p - no) ( cTp-£- j |( 5 np|^ = C -I- ( 5 up 5 nn) 


(34) 


(35) 


(36) 


where we have introduced the pattern speed (Tp. From this expression we see that the imaginary part of the left-hand 
side must vanish, so we should have Im C = Q. Allowing the pattern speed to be complex, i.e. using (jp = an -\- iaj we 
find that the following condition must be satisfied 


Im C — anai 


2 - 


vq 

(Tfl 


1 - 


(Tp — Ho 




+ 


2 - 


no 

(Tfl 


1 + 


\5vp 


= 0 . 


(37) 


If we are to have an unstable mode, aj yf 0, the frequency clearly must be such that the factors multiplying the absolute 
values of the two velocities have different signs. 

Let us first consider the case when the factor multiplying |5nn|^ is negative. Then we find that an instability is only 
possible if (Tij/uo < 0, and the following condition is satisfied: 
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Figure 2. An illustration of the instability criterion ( |40| ) which is relevant for the example considered in Section IIB. This 
example is constructed by introducing 2 = ct^/uq, and then showing the four curves: 0 ( 2 ) = 2 , b{z) = (1 + 7 ^/ 2 ^)/ 2 , c{z) = 1 and 
d(z) = (2 — 1)^. Here we have taken c^/vq = 1 and 7 ^ = Cp/c^ = 0.0379. Criterion (kol) is satisfied when d < c and a < b (in the 
grey area). The corresponding range is well approximated by 0 < o-^/vq <1/2. 


0 < 


<yR 

Vo 


1 

<2 


Fp - Wol 


-1 


This shows that we must have 
Cn 

I Wo 


> 



> 

^-1 

_ 

(7R 

Vo 


Wo 


Vo 


+ 1 > 1 


(38) 


(39) 


which constrains the permissible frequencies to the range |CTi?| < Cn — |wo|- Thus we see that the flow must be subsonic, 
i.e. |wo| < Cn. 

In the case when the factor multiplying |5wp|^ is negative we can only have an instability if an/vo > 0. We also 
require 


«<-<! 


or 

1 

2 


Wo 


1 - 


1 + 


< 


1 + ■ 


CTfl - Wo 


1 - 


Up - Wo 


if jup - 

vof < cl 

UR 1 


h 

V 

V 


Vo 2 

V 


if 


Up - Wo > Cn 


(40) 


(41) 


For the example illustrated in Figure the condition that must be satisfied is (^0|). It is useful to notice two things 
about this criterion. First of all, any unstable mode for which an > Cp must lie in the range 0 < Cp < < wq. Secondly, 

when (T_R >> Cp the permissible range will be well approximated by 0 < au < wo/2, cf. Figure |^. As is clear from 
Figure ^ the unstable modes satisfy this last, and most severe, criterion. 

It is worth noting that the instability can be discussed in terms of a simple energy argument [see |Casti etlii (1998) 
and Pierce (1974) for similar discussions in other contexts]. After averaging over several wavelengths, the kinetic energy 
of the protons is 


Sp 




Meanwhile we get for the neutrons 
En ^ ^(rin + 6nn){vo + (5wn)^ 


which leads to 
m^rin 


En: 


2 LO + kvo r 2 

Wo H-- dVn 

LO — KVo 


(42) 


(43) 


(44) 


from which we see that the energy in the perturbed flow is smaller than the energy in the unperturbed case, which 
means that we can associate the wave with a “negative energy”, when 


© 0000 RAS, MNRAS 000, 000-000 
































The superfluid two-stream instability and pulsar glitches 7 


a 



0.2 


0.4 


0.6 


0.8 


Figure 3. An illustration of the various domains of instability for the simple “analytic” equation of state of Prix et al (200!^). An 
absolute instability (see discussion in the main text) is active in the grey areas (also labeled I and II). The two-stream instability 
is, in principle, relevant in the remaining parameter space. The dashed curves indicate the onset of instability when the relative 
flow is equal to the neutron sound speed (y = 1). For slower flows, these critical curves approach the absolute instability regions. 
The region where the two-stream instability may operate in physical flows therefore lies between each dashed curve and the nearest 
grey area. For comparison we also indicate the curve in the Xp — <j plane traced out by the PAL equation of state (discussed in 
Section HE) as the density is varied from that near the crust/core interface (u = 0.64) to flve times that of nuclear saturation 
(u = 5). 


— vo<-^<vo, i.e. — Vo < (JR < Vo ■ (45) 

A wave that satisfies 0 < au < vq moves forwards with respect to the protons but backwards according to an observer 
riding along with the unperturbed neutron flow. As we have seen above, the unstable modes in our problem satisfy this 
criterion and hence it is easy to explain the physical conditions required for the two-stream instability to be present. 


2.4 Results for a simple model equation of state 


Having established that the two-stream instability may be present in superfluids, we want to assess to what extent one 
should expect this mechanism to play a role for astrophysical neutron stars. To do this we will consider two particular 
equations of state. The results we obtain illustrate different facets of what we expect to be a rich problem. 

We begin by making contact with our recent analysis of rotating superfiuid models (Prix et al 2002) as well as the 
study of oscillating non-rotating stars by Prix & Rieutord (2002). From the definitions above we have 


p ‘-^np 
■n 


(46) 

(47) 


We combine these results with the explicit structure matrix given in eq. (144) of Prix et al (2002), which is based on a 
simple “analytic” equation of state. This leads to 

2 n-D 


Ip 2 
-(T 


= 


1-p 


cr(l — 2 Xp) 

1 Xp 


where Xp = np/(np -P rip) is the proton fraction, and cr is defined by 


(48) 

(49) 

(50) 


As discussed by Prix et al (2002), o is related to the “symmetry energy” of the equation of state, cf. Prakash et al 
(|198^). 
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8 N. Andersson, G.L. Comer and R. Prix 


The instability regions for this model equation of state are illustrated in Figure]^. A key feature of this figure is the 
presence of regions of “absolute instability”. This happens when a? > . Then there exist unstable solutions already 

for vanishing background flow, y = 0. That this is the case is easy to see. Consider in the limit 3 / = 0. In the limit 
we can solve directly for x^\ 


X 


2 



(51) 


from which it is easy to see that one of the roots for x'^ will be negative if a? > b^. Hence, one of the roots to the quartic 
dd^ ) will be purely imaginary. 

The physics of this instability is quite different from the two-stream instability that is the main focus of this paper. 
Yet it is an interesting phenomenon. From the above relations we find that corresponds to 




(52) 


In the discussion by Prix et al ( P002|) it was assumed that “reasonable” equations of state ought to satisfy this condition. 
We expected this to be the case since the structure matrix would not be invertible if its determinant were to vanish at 
some point. We now see that this constraint has a strong physical motivation: The condition is violated when > b^, 
i.e. when we have an absolute instability. The regions where this instability is active are indicated by the grey areas in 
Figure H. 


2.5 Results for the PAL equation of state 


In order to strengthen the argument that the two-stream instability may operate in astrophysical neutron stars, we have 
considered a “realistic” equation of state due to Prakash, Ainsworth and Lattimer (PAL) (1988). The advantage of this 
model is that it is relatively simple. In particular, it leads to analytical expressions for the various quantities needed in 
our analysis. The energy density of the baryons for the PAL equation of state can be written 

S{nn, Up) = (rin + rip) [Eo{u) + S{u){l - 2xp)'^] , (53) 

where Eq corresponds to the energy per nucleon, S corresponds to the “symmetry energy” (and is closely related to a 
in the “analytic” equation of state discussed above), and u = (rip -|-np)/no with no = 0.16 fm“® the nuclear saturation 
density. Eq takes the following form: 

2 

Eo{u) = -I- Bqu -I- CqvE -I- 3 Cia~^ — arctan , (54) 

i=l 


with Ao = 22.11 MeV, Bo = 220.47 MeV, Co = -213.41 MeV, cr = 0.927, Ci = -83.84 MeV, C2 = 23.0 MeV, ai = 2/3, 
and 02 = 1/3. The symmetry term is 

S{u) = As [u"/® - P(m)] -b SoE{u) , (55) 


with As = 12.99 MeV and So = 30 MeV. Here F{u) is a function satisfying F{1) = 1 which is supposed to simulate 
the behaviour of the potentials used in theoretical calculations. In our study we have only considered F{u) = u, which 
is one of four possibilities discussed by Prakash et al (1988). 

We further need to account for the energy contribution of the electrons, which is important since the electrons are 
highly relativistic inside neutron stars. Hence, they can obtain high (local) Fermi energies which may be comparable 
with the proton (local) Fermi energies. Considering only the electrons, the leptonic contribution to the energy density 
is given by (in units where the speed of light is unity) (Shapiro & Teukolsky 1983) 


me / F\ 

C — ^3 X (Xe ) 


(56) 


where me = mb/1836 is the electron mass (in terms of the baryon mass mu), Ae 
wavelength, and 


Xix) = (^x [1 -b 2;^] [1 -b 2 a;^] — In a: -b (l -b x^) ^ 


hi me is the electron Compton 


(57) 
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Figure 4. Two-stream instability results for the PAL equation of state. Left panel: The region where the two-stream instability 
is present (grey area) is shown as a function of the density parameter u. We indicate the location of the core-crust boundary 
{u 0.64) by a vertical dashed line. Our model is only relevant for the core fluid, i.e. to the right of the vertical line. Finally, the 
horizontal dashed line indicates when the relative flow is equal to the (neutron) sound speed. We expect the superfluid degeneracy 
to be broken beyond this level, so an instability located above this line is unlikely to have physical relevance. The results indicate 
that there may be a region of instability immediately below the crust. Right panel: The corresponding oscillation frequencies. 
Particularly notable is the point near u = 2 where the two critical curves cross. At this point the symmetry energy cr changes sign, 
cf. Fig. Id and there exist a particular density such that the two fluids are uncoupled, cf. (pl[). 


Xe = 1836 




1/3 


1/3 


(58) 


In doing this calculation we have assumed local charge neutrality, i.e. rie = rip. The above energy term is added linearly 
in the equation of state. 

Having obtained an expression for the total energy as a function of the density, we can derive explicit expressions 
for all quantities needed to discuss the two-stream instability. First we need to determine the proton fraction Xp. We do 
this by assuming that the star is in chemical equilibrium, i.e. 


fT = fT ■ (59) 

Solving ( ^ ) for Xp provides us with the proton fraction as a function of u. Given this, and the relevant partial derivatives 
oi £ -\- Eg we can readily evaluate the symmetry energy and well as the sound speeds cf, Cp and the chemical coupling 
parameter C. With this data we can determine the two parameters and which are needed if we want to solve the 
local dispersion relation (^^. The results we obtain for the proton fraction and the symmetry energy are indicated in 
Figure ^ We consider the range from u = 0.64, presumed to correspond to the core-crust boundary, to u = 5 which 
represents the deep core of a realistic neutron star. The corresponding results for the two-stream instability are shown 
in Figure ^ From this hgure we can see that the two-stream instability may operate (albeit at comparatively large 
relative flows) in the region immediately below the crust. Finally, we hnd that the conditions at the core-crust interface 
are such that = 0.0249 and = 0.0379. These are the values we chose for the example in Section IIB and hence the 
results shown in Figure ^ correspond to a physically realistic model. 


3 ROTATING SHELLS: A CONNECTION WITH PULSAR GLITCHES? 

The fact that a superfiuid neutron star may, in principle, exhibit the two-stream instability does not necessarily prove that 
this mechanism will be astrophysically relevant. Yet, it is an intriguing possibility given that the mechanism underlying 
the enigmatic glitches observed in dozens of pulsars remains poorly understood. One plausible astrophysical role for the 
superfluid two-stream instability would be in this context: Perhaps this instability serves as trigger mechanism for large 
pulsar glitches? 
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3.1 Dispersion relation for rotating shells 

Our aim in this Section is to construct a toy problem that allows us to investigate a possible connection between the 
two-stream instability and pulsar glitches. A suitably simple problem corresponds to two fluids, allowed to rotate at 
different rates, confined within an infinitesimally thin spherical shell. By assuming that the shell is infinitesimal we 
ignore radial motion, i.e. we restrict the permissible perturbations of this system is such a way that the perturbed 
velocities must take the form 

1 , .S ^ rm. ^ ^ 1 + lW. ir 


Svx = - 


p. ,UUt)d^Yreg + ^UL{t)deYr 

JX sm u IX 


(60) 


where {6, ip) are the spherical harmonics and R is the radius R of the shell. This means that the system permits 
only toroidal mode-solutions. In other words, all oscillation modes of this shell model are closely related to the inertial 


r-modes of rotating single fluid objects (Papaloizou & Pringle 


1978 


Provost et al 1981). 


The perturbation equations for the conhguration we consider have been derived in a different context and the 


complete calculation will be presented elsewhere (Andersson et al 2002). Our primary interest here concerns whether 
the modes of this system may suffer the two-stream instability. The presence of the instability in this toy problem would 
be a strong indication that it will also be relevant when the shells are “thick” and radial motion is possible. That is, 
when we consider a rotating star that contains a partially decoupled superfluid either in the inner crust or the core. 

As discussed in Section IIC, the two-stream instability can be understood in terms of negative energy waves. In the 
current problem, the criterion for waves to carry negative energy according to one fluid but positive energy according 
to the other fluid is that the mode pattern speed [we are assuming a decomposition expiiuit -\- irmp)] 

LO 

Cp = 


(61) 


lies between the two (uniform) rotation rates. In other words, a necessary condition for instability is 


(62) 


where we have assumed that the superfluid neutrons lag behind the charged component, as is expected in a spinning 

down pulsar. _ _ 

After a somewhat laborious calculation, see Andersson et al ( |2002 ) for details, one finds that the dispersion relation 
for the toroidal two-fluid modes of the shell problem is 


{ — I (I -t- 1) (1 — £n) “t mfln] + mSn [2 — / (/ -t- 1)] [rip — fin]} 

X {—/ (I + 1) (1 — Ep) [lo -|- mQp] + 2mQp -1- mep [2 — I {I + 1)] [fin — flpj} 
— {I {I + 1))^ EnEp (a; -I- mflu) [oJ + mflp} = 0 . 


(63) 


This equation should be valid under the conditions in the outer core of a mature neutron star, where superfluid neutrons 
are permeated by superconducting protons. 

We rewrite this dispersion relation in terms of the entrainment parameter used by Prix et ^ (2002), namely 
Ep = ^ , (64) 


e = Ep 


Pp 


the frequency as measured with respect to the rotation of the protons, 
UJ 


K = 


mflp 


and a dimensionless measure of the relative rotation, 
fin 


A = 


flp 


(65) 


( 66 ) 


With these definitions we get 

{/(I -t- 1)[1 — ^p(l T e')](k -f A) — 2(1 — Xp)A -t- Xpsil — !)(/ -t- 2)(1 — A)} 

X {l{l + 1)[1 - s]{k +l)-2-e{l-l)il + 2)(1 - A)} 

— [l{l + l}]^Xpe^(K -I- A)(k -I- 1) = 0 . 

In terms of these new variables a mode would satisfy the necessary instability criterion (^ 
— A < K < —1 . 


if K is such that 


(67) 


( 68 ) 
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Figure 5. Instability regions for the shell problem. We show results for 1 = 2 and relative rotation such that the neutrons rotate 
90% slower (left frame) or faster (right frame) than the protons. The two-stream instability operates in regions 1-2. The grey box 
corresponds to the “physically reasonable” part of parameter space (for the core of a neutron star). 


As we will now establish, there exist modes that satisfy this criterion for reasonable parameter values. 


3.2 An extreme example: The quadrupole modes 

Let us first consider the case of quadrupole oscillations, i.e. take I = 2. Typical results for this case are shown in Figures ^ 
and ^ The first figure illustrates the regions of the Xp — e parameter space for which an instability is present in the case 
when (i) the neutrons rotate at a rate that is 90% faster than that of the charged component, and (ii) the neutrons lag 
behind by the same fraction. The second figure shows the mode-frequencies corresponding to the second case [simply 
obtained by solving the quadratic for the particular value Xp = 0.05. This figure shows the presence of unstable 
modes within the range of values for e that we take to be physically realistic (Prix et al 2002): 0.4 < e < 0.7. From 
this figure we immediately deduce two things. First, we see that the unstable modes indeed satisfy the criterion (|68|). 
Secondly, the unstable modes may have imaginary parts as large as Im k ~ 0.15. It is useful to ask what this implies 
for the growth time of the instability. Since our results only depend on the azimuthal index m through the scaling 


Im uj = mflp Im k 


(69) 


we see that the fastest growth time corresponds to the m = I modes. The e-folding time for the instability follows from 
1 1 


te = 


27rm Im k 


P s 


(70) 


where P represents the observed rotation period of the pulsar (presumably corresponding to the rotation of the charged 
component, i.e. P/2ti/Q, p). 

The above example shows that the two-stream instability does indeed operate in this shell problem. In fact, the 
analysis goes beyond the local analysis of the plane parallel problem in Section IIB since we have now solved for 
the actual unstable modes (satisfying the relevant boundary conditions). We see that, as expected for a dynamical 
instability, the growth time of an unstable mode may be very short. However, the relative rotation rates required to 
make the quadrupole modes unstable in the range 0 < Xp < 1 and 0 < e < 1 are likely far too large to be physically 
attainable. In this sense the results shown in Figs. are, despite being instructive, somewhat extreme. 


3.3 A realistic mechanism for triggering pulsar glitches? 

A quantity of key importance for this discussion is the rotational lag between the two components. In order to be able 
to argue that the two-stream instability is relevant for pulsar glitches we need to consider lags that may actually occur 
in astrophysical neutron stars. To estimate the size of the rotational lag required to “explain” the observed glitches we 
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Figure 6. Real and imaginary parts (left and right frame, respectively) of the mode frequencies for / = 2, A = 1.9 and Xp = 0.05. 
The corresponding instability region can be deduced from the right panel of Fig. H. The grey area in the left frame indicates the 
region where an instability is permissible according to (|68|). 


assume that a glitch corresponds to a transfer of angular momentum from a partially decoupled superfiuid component 
(index n) to the bulk of the star (index p). Then we have 

/nlAfinl « /pA^p ^ A^p « :^|Ann| (71) 

where Ix are the two moments of inertia. Now assume that the decoupled component corresponds 1% of the total 
moment of inertia, eg. the superfluid neutrons in the inner crust or a corresponding amount of fluid in the core. This 
would mean that In ~ 10“^/p, and we have 

Afip « lO'^IAfinl . (72) 


Combine this with the observations of large Vela glitches to get 


Aflp 


10 ' 


lAfln 


10 “ 


flp flp 

In other words, we must have 

|Afln| « 10"‘‘np . 


(73) 


(74) 


If we assume that the glitch brings the two fluids back into co-rotation, then we have Afln = fin — flp = All and we see 
that the two rotation rates will maximally differ by one part in 10"^ or so. Rotational lags of this order of magnitude have 
often been discussed in the context of glitches. Even though the key quantity in models invoking catastrophic vortex 
unpinning in the inner crust — the pinning strength — is very uncertain, and there have been suggestions that the 


pinning force is too weak to allow a build up of the required rotational lag (Jones 1998), typical values considered are 
consistent with our rough estimate. In addition, frictional heating due to a difference in the rotation rates of the bulk 
of a neutron star and a superfluid component has been discussed as a possible expla nation for the fact that o l d isolated 
pulsa rs seem to be somewhat hot ter than expected from standard cooling models ( Shibazaki fc Lamlj 1989 ; Larson fc 
Lin^ 199£ ). Larson fc Link ( 1991 ) argue that a lag of 


AQ 

aT 


(3.2 X 10~'‘ - 9.5 X 10“®) X (-^) 


(75) 


could explain the observational data. Finally, the presence of rotat ional lags of the proposed magnitude is supported by 
a statistical analysis of 48 glitches in 18 pulsars ( Lyne etii] 2000). This study suggests that the critical rotational lag 
at which a glitch occurs is 


All 

flp 


5 X 10“ 


(76) 
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Figure 7. The two-stream instability regions (labeled 1 and 2) for the case when the superfluid neutrons lag behind the 
superconducting protons in such a way that Af2/f7p = 5 X 10“^ or A = 1.0005. We show results for four different values of 1. 
There are no unstable modes in this part of parameter space unless I > 65. Recall that Xp < 0.1 and 0.4 < £ < 0.7 would be 
reasonable parameter values for a neutron star core. The grey diagonal line represents the singularity discussed in Section HID 
and the dashed curves indicate the instability regions in the large I limit. 


In order to make our shell model problem more realistic we consider the case when the superfluid neutrons lag 
behind the superconducting protons, and take Afl/flp = 5 x 10“^, or A = 1.0005, as a representative value. With this 
rotational lag, we see from dosl ) that the unstable modes must be such that —1.0005 < k < —1. A series of results for 
this choice of parameters are shown in Figures 

Figure ^ illustrates the fact that, if we decrease the rotational lag then the two-stream instability will not be active 
(in the interesting region of parameter space) for low values of 1. For a smaller rotational lag the instability acts on a 
shorter length scale. For A = 1.0005 we And that we must have Z > 65 in order for there to be a region of instability in 
the part of the Xp — e plane shown in Figure This means that the instability only operates on length scales shorter 
than -kR/I « 500 m (if we take the shell radius to be i? = 10 km). 

Figure ^ shows the real part of the mode frequencies for I = 100 and various values of Xp. From this figure we can 
see that the instability always occur in the region suggested by the instability condition (^), i.e. two real frequency 
modes never merge to give rise to a complex conjugate pair of solutions outside the grey areas indicated in the various 
panels of Figure The imaginary parts for I = 100 and the various values of Xp considered in Figure ^ are shown 
in Figure p. From this figure we see that the imaginary part of k typically reaches values of order 10~'*. In fact, by 
comparing similar results for different values of I we have found that the largest attainable imaginary part of k varies by 
less than one order of magnitude as I increases from 100 to 1000. Thus we estimate that the typical instability growth 
time for A = 1.0005 will be 


te 


IQ-* P 
m 2 -k 


(77) 


For a star rotating at the rate of the Vela pulsar, P = 89 ms, we would have te « 1.4 s for I = m — 100. Interestingly, 


this predicted growth time is significantly shorter than the resolved rise time of a large Vela glitch tgutch < 40 s (Dodson 
et ji| ^002| ). 


3.4 Approximate results: The large I limit 

The results obtained above indicate that the two-stream instability is likely to act on modes with rather short wave¬ 
lengths. Given this it makes sense to consider the large I limit in more detail. Keeping only the leading order term 
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Figure 8. Mode frequencies (represented by Re k.) as function of the entrainment parameter e for I = 100, A = 1.0005 and various 
values of Xp. The range in which an instability is permissible [according to (S)] is indicated by the grey areas. 



Figure 9. The imaginary part Im k. (determining the growth rate) of the unstable modes for A = 1.0005 and I = 100 is shown 
as a function of e for several Afferent values of Xp (as indicated in the figure). Essentially, these results represent horizontal cuts 
through the frames in Figure M. 


(proportional to /^) in (jo^ we have the dispersion relation 

(1 — £ — Xp)f^^ + [1 — 2s — Xp + 2xpe + (1 — a;p — 2xpe)A] k, 

+ {1 — 2e — Xp)A + eXp + £(1 — Xp)A^ — 2xpS^{l — A)^ = 0 . 

We can readily write down the solutions to this quadratic: 

= 7 - 7 :;- - -r {1 — 2s — Xp 2exp + (1 — Xp — 2 xp£)A 

2(1 — Xp — s) 

± (1-A)[(l-2e)(l-2e-2;p)(l-a:p-2£a;p)]^/q . 


(78) 


(79) 
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It turns out that we can learn a lot about the problem from this expression. The most obvious feature is the fact that 
K± have a singularity when 1 — a;p — e = 0. It is straightforward to show that one of the roots will become infinite at 
this point, while the second root becomes: 

1 


l-2x. 


[A — Xp(l + A) + 2a;pe(l — A)] . 


(80) 


For the examples shown in Figure Q this special case corresponds to the intersection between the instability domain and 
the diagonal line 1 — Xp — e = 0, as indicated in Figure |^. 

It is also straightforward to deduce the regions of instability from the sign of the argument of the square-root. We 
see that we will have an unstable mode when 
1 . , „ 1 


e < 


and 


1 — 2£ < Xp < 


l + 2e 


Region 1 


(81) 


£ > 


1 


and 




Region 2 . 


(82) 


2 ^ 1 -\- 2 £ 

These regions are also indicated in Figure Out of these two possible instability regions, the first is most likely to be 
relevant for neutron stars since it allows the instability to be present already for small proton fractions. 

Finally, we can use (^^ to show the existence of an extremum at 

iVrf + 7 


£ = i(^-l)^0.39 


4^17 + 3 


0.39 . 


The corresponding imaginary part of k would represent the fastest possible growth time for an unstable mode located 
in region 1 of Figure This leads to the estimate 


Im K < 0.24(A - 1) 

or an estimate of the shortest growth time: 
t : 


6.7 X 10"^ , 


I'V ^ 1 

1 1 

Up J 

VO.l s) 


(83) 


(84) 


which agrees well with the result we previously obtained for I = m = 100. Hence, this simple formula can be used to 
estimate the fastest growth rate of the two-stream instability in our shell model for different parameter values. 


4 DISCUSSION 


In this paper we have in troduced the superflu id two-stream instability: A dynamical instability analogous to that known 
to operate in plasmas ( Anderson et al |2001 ), which sets in once the relative flow between the two components of the 
system reaches a critical level. We have studied this instability for two different model problems. First we analysed a 
local dispersion relation derived for the case of a background such that one fluid was at rest while the other had a 
constant flow rate. This provided a proof of principle of the existence of the two-stream instability for superfluids. Our 
analysis was based on the two-fluid equations that have been used to model the dynamics of the outer core of a neutron 
star, where superfluid neutrons are expected to coexist with superconducting protons and relativistic electrons. These 
equations are analogous to the Landau model for superfluid Helium and should also (after some modifications to 
incorporate elasticity and possible vortex pinning) be relevant for the conditions in the inner crust of a mature neutron 
star. Thus we expect the two-stream instability to be generic in dynamical superfluids, possibly limiting the relative 
flow rates of any multi-fluid system. Our second model problem concerned two fluids confined within an infinitesimally 
thin spherical shell. The aim of this model was to assess whether the two-stream instability may be relevant (perhaps 
as trigger mechanism) for pulsar glitches. The results for this problem demonstrated that the entrainment effect could 
provide a sufficiently strong coupling for the instability to set in at a relative flow small enough to be astrophysically 
plausible. Incidentally, the modes that become dynamically unstable in this problem are the superfluid analogues of 


^ Even though we do not discuss this issue in detail here, it is exciting to contemplate possible experimental verification of the 
superfiuid two-stream instability in, for example, superfluid “^He. 
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the inertial r-modes of a rotating single flnid star. This is interesting since the r-modes are known to be secnlarly 
unstable due to the emission of gravitational radiation (Andersson & Kokkotas 2001). In fact, the connection between 
the two instabilities goes even deeper than this since the radiation-driven secular instability is also a variation of the 
Kelvin-Helmholtz instability. In that case, the two fluids are the stellar fluid and the radiation. 

In order for an instability to be relevant the unstable mode must grow faster than all dissipation timescales. In the 
case of a superfluid neutron star core the main dissipation mechanisms are likely to be mutual friction and “standard” 
shear viscosity due to electron-electron scattering. Since we have tried to build a plausible case for the two-stream 
instability to be relevant for pulsar glitches we would like to obtain some rough estimates of the associated dissipation 
timescales. To do this we first use an estimate of when mutual friction is likely to dominate the shear viscosity [due to 
Mendel^ (1991b)]: 

2 


n > 100 


10 ® cm 
A 




where we assume that the wavelength of the mode is A = ttR/I. We can write this as 


P < 0.62r 


f——1 

Vio^ kJ 


(85) 


( 86 ) 


which suggests that the shear viscosity will be the dominant dissipation mechanism for large values of 1. For example, 
for a neutron star rotating with the period of the Vela pulsar mutual friction would dominate for / < 15 or so (assuming 
Tt ~ 5). To estimate the shear viscosity damping we can use results obtained for the secular r-mode instability. In 
particular, Kokkotas & Stergioulat (1999) have shown that for a uniform density star one has 


M 


3 M 


4Tr{l-l){2l + 3)'nR SnPriR’ 

If we use the shear viscosity coefficient for electron-electron scattering 




6 X 10" 


10 ^® g/cm'" 


10® K 
T 


, .nisf M V/lO km\® /lO® kV , 

g/cm s-2.7x10 g/cii 


we get 


1.2 X 10 "^ 


/ 1-4M0 \ / P / r y 
V M AlOkm/ UO^K/ 


(87) 


( 88 ) 


(89) 


P \ M /VlOkm/ VlO^Ky 
We want to compare this damping timescale to the growth rate of the unstable modes in our shell toy-model. 
Combining ( ^^ with dsi] ) we estimate that in order to have an instability we must have 

Afl/rip ^ /n 1 / 1 ^ V / D x5/ rr x2 


I < 90 


5 X 10-4 


( 0-1 s\ / 1.4M0 \ ( R y ( T y 

\ P )\ M AlOkmy* UO'^ Ky* ■ 


(90) 


Let us now consider the case of the Vela pulsar. Estimati ng the core temperature as 5 x 10^ K [roughly two orders of 
magnitude higher than the observed surface temperature (|Page| |l997 )[ we deduce that only modes with I > 2500 or so 
are likely to be stabilized by shear viscosity. Given that our results indicate that the two-stream instability is active 
for much smaller values of I, cf. the results shown in Figure ^ we conclude the dissipation is unlikely to suppress the 
instability in sufficiently young neutron stars. Incidentally, the length scale corresponding to a mode with I — 2500 
would be about ten meters. This is an interesting result since one can show t hat a large glitch could be explained by a 
small fraction (~ 10”“^) of the neutron vortices moving a few tens of metres ( Cordes et al 198j ). 

Obviously, the situation changes as the star cools further. Based on the above estimates one can show that shear 
viscosity will suppress all modes with I > 65 (i.e. all unstable modes for the case considered in Fig. ^) if the core 
temperature is below 8 x 10® K. This means that the two-stream instability may not be able to overcome the viscous 
damping in a sufficiently cold neutron star, which is consistent with the absence of glitches in mature pulsars. 

We believe that the results of this paper suggest that the superfluid two-stream instability may be relevant in the 
context of pulsar glitches. If this is, indeed, the case then what is its exact role? The answer to this question obviously 
requires much further work, but it is nevertheless interesting to speculate about some possibilities. Most standard models 


for glitches are based on the idea of catastrophic vortex unpinning in the inner crust (Anderson & Itoh 1975). This 
is an attractive idea since the glitch relaxation (on a timescale of days to months) would seem to be well described 
by vortex creep models (Cheng et al 1988). An interesting scenario is provided by the thermally induced glitch model 
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discussed by Link & Epstein (1996). They have shown that the deposit of 10“^^ erg of heat would be sufficient to induce 
a Vela type glitch. The mechanism that leads to the unpinning of vortices, eg. by the deposit of heat in the crust, is 
however not identified. We believe that the two-stream instability may fill this gap in the current theory. It should, of 
course, be pointed out that glitches need not originate in the inner crust. In particular, Jones (1998) has argued that 
the vortex pinning is too weak to explain the recurrent Vela glitches. If this argument is correct then the glitches must 
be due to some mechanism operating in the core fluid. Since the model problems we have considered would be relevant 
for the conditions expected to prevail in the outer core of a mature neutron star, our results show that the two-stream 
instability may serve as a trigger for glitches originating there. The key requirement for the instability to operate is the 
presence of a rotational lag. It is worth pointing out that such a lag will build up both when there is a strong coupling 
between the two fluids (i.e. when the vortices are pinned) and when this coupling is weak. One would generally expect 
the strength of this coupling to vary considerably at various depths in the star (Langlois et al 1998), and it is not yet 
clear to what extent a rotational lag can build up in various regions. This is, of course, a key issue for future theoretical 
work on pulsar glitches. 

One final relevant point concerns the recent observation of a Vela size glitch in the anomalous X-ray pulsar IRXS 
J170849.0-400910 (Kaspi et a] 2000| ). This object has a spin period of 11 s, which means that any feasible glitch model 
must not rely on the star being rapidly rotating. What does this mean for our proposal that the two-stream instability 
may induce a glitch? Let us assume that the rotational lag builds up at the same rate as the electromagnetic spindown 
of the main part of the star (i.e. that the superfluid component does not change its spin rate at all under normal 
circumstances). Then the lag would be Af2 « tO. after time t. If there is a critical value at which a glitch will happen 
(corresponding to Aflcrit) then the interglitch time tg could be approximated by 

AUcrit U n Aflcrit 


where r is the standard “pulsar age”. This argument implies the following: i) For Aflcrit/U « 5 x 10”'* we would get 
tg « 10”®r. This (roughly) means that only pulsars younger than 10"^ yr would be seen to glitch during 30 years of 
observation, which accords well with the fact that only young pulsars are active in this sense, ii) There is no restriction 
on the rotation rate in this scenario; a star spinning slowly may well exhibit a glitch as long as its spindown rate is fast 
enough. This means that one should not be surprised to find glitches in stars with extreme magnetic fields (magnetars). 

This paper is only a first probe into what promises to be a rich problem area. Future studies must address issues 
concerning the effects of different dissipation mechanisms, the nonlinear evolution of the instability, possible experimental 
verification for superfluid Helium etcetera. These are all very interesting problems which we hope to investigate in the 
near future. 
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